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ABSTRACT 

When a star is tidally disrupted by a supermassive black hole (BH), roughly half of its mass 
falls back to the BH at super-Eddington rates. Being tenuously gravitationally bound and un¬ 
able to cool radiatively, only a small fraction /in ^ 1 of the returning debris will likely be 
incorporated into the disk and accrete, with the vast majority instead becoming unbound in an 
outflow of velocity ~ 10^ km s“^ This slow outflow spreads laterally, encasing the BH. For 
months or longer, the outflow remains sufficiently neutral to block hard EUV and X-ray radia¬ 
tion from the hot inner disk, which instead becomes trapped in a radiation-dominated nebula. 
Ionizing nebular radiation heats the inner edge of the ejecta to temperatures of T > few 10"^ K, 
converting the emission to optical/near-UV wavelengths where photons more readily escape 
due to the lower opacity. This can explain the unexpectedly low and temporally constant ef¬ 
fective temperatures of optically-discovered TDE flares. For BHs with relatively high masses 
M, > 10^Mo the ejecta can become ionized at an earlier stage, or for a wider range of view¬ 
ing angles, producing a TDE flare which is instead dominated by thermal X-ray emission. 
We predict total radiated energies consistent with those of observed TDE flares, and ejecta 
velocities that agree with the measured emission line widths. The peak optical luminosity for 
M, < 10 ^Mq is suppressed as the result of adiabatic losses in the inner disk wind, possibly 
contributing to the unexpected dearth of optical TDEs in galaxies with low mass BHs. In the 
classical picture, where /n ~ 1, TDEs de-spin supermassive BHs and cap their maximum 
spins well below theoretical accretion physics limits. This cap is greatly relaxed in our model, 
and existing Fe Ka spin measurements provide suggestive preliminary evidence that /n < 1. 


1 INTRODUCTION 

Occasionally a star in a galactic nucleus is perturbed onto a highly 
eccentric orbit that brings it close enough to the central supermas¬ 
sive black hole (BH) to be torn apart by tidal forees (Hills 1975; 
Carter & Luminet 1983). The rapid accretion of stellar debris fol¬ 
lowing such a tidal disruption event (TDE) was predicted to power a 
luminous electromagnetic flare, detectable to large distances (Rees 
1988). 

TDE flares have been discovered at hard X-ray (e.g.. Bloom 
et al. 2011, Burrows et al. 2011; Cenko et al. 2012b, Pasham et al. 
2015), soft X-ray (e.g.. Bade et al. 1996, Grupe et al. 1999, Ko- 
mossa & Bade 1999, Komossa & Greiner 1999, Donley et al. 2002, 
Esquej et al. 2008, Maksym et al. 2010, Saxton et al. 2012), ultravi¬ 
olet (e.g.. Stern et al. 2004, Gezari et al. 2006, 2008, 2009), and op¬ 
tical (van Velzen et al. 2011, Gezari et al. 2012, Cenko et al. 2012a, 
Arcavi et al. 2014, Chornock et al. 2014, Holoien et al. 2014, Vinko 
et al. 2015) wavelengths (see Komossa 2015 for a recent review). 
These events are generally characterized by angular positions con¬ 
sistent with the centers of their host galaxies; luminosities close to 
or exceeding the Eddington luminosity for the inferred BH mass; 
and a characteristic light curve decay similar to the predicted oc 
fail-back rate of stellar debris following a TDE (Rees 1988; Ayal 
et al. 2000; Lodato et al. 2009). 

Despite some similarities, it is not yet clear how this heteroge¬ 
neous class of events, detected by surveys at different wavelengths 
with different selection criteria and biases, fits together into a co¬ 


herent picture of the TDE phenomenon. For instanee, while some 
jetted TDEs (those characterized by nonthermal X-ray and syn¬ 
chrotron radio emission) show luminous thermal optical emission 
(Pasham et al. 2015), most optical and soft X-ray discovered flares 
are radio quiet (Bower et al. 2013; van Velzen et al. 2013). The first 
TDEs were detected in soft X-rays, yet many optically-discovered 
events produce no detectable X-ray emission (e.g. PSl-lOjh, Gezari 
et al. 2012, and ASAS-SN14ae, Holoien et al. 2014), or do so only 
after a significant delay (e.g., D3-13, Gezari et al. 2006 and Dl-9, 
Gezari et al. 2009). 

The peak effective temperature of thermal radiation pro¬ 
duced by radiatively-efficient accretion onto a BH of mass M, = 
lO^M. 6 Mo is given by (e.g., Miller 2015) 
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where M is the accretion rate normalized to the Eddington rate of 
Mfidd = AnGM.I k^s ~ 0.4 cm^ g"^ is the electron scat¬ 
tering opaeity, rj = Lacc/Mc^ is the radiative effieiency, Lace is the 
accretion luminosity, and Rin = 6Rg is the radius of the innermost 
stable orbit for a non-rotating BH in units of the gravitational radius 


Rg = GM./c/ 


One of the puzzles of the optically-discovered TDEs is why 
their measured color temperatures are consistently much lower than 
those predicted by standard aceretion models. The best-fit Planck 
temperatures are usually in the range of^ 2-5x10'^K (van 
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Velzen et al. 2011, Gezari et al. 2012, Arcavi et al. 2014, Chornock 
et al. 2014, Holoien et al. 2014, Vinko et al. 2015), well below the 
predicted values of Teff ^ 3xl0^-10^K from equation (1) for 
M, < 10^ K and the range of accretion rates, M > MEdd. thought 
to characterize TDEs. PSl-lOjh maintained a constant optical-UV 
spectral slope, even as the total flux decreased by more than a factor 
of 10 (Gezari et al. 2012). 

The low temperatures of TDE flares have generally been at¬ 
tributed to the presence of a “reprocessing layer,” which absorbs 
higher energy radiation from the inner disk and reradiates it at 
lower frequencies (Loeb & Ulmer 1997, Guillochon & Ramirez- 
Ruiz 2013). If such a layer is in thermal equilibrium, then to ex¬ 
plain the observed temperatures the layer should be located at radii 
of > 10^ - lO^Rg, which are much larger than the radius < lOORg 
at which the stellar debris will circularize into the accretion disk. 

One possible source of radially-extended matter are collisions 
between the eccentric streams of stellar debris created during the 
disruption, which may contribute to the process of debris circu¬ 
larization (e.g., Hayasaki et al. 2013, Cheng & Bogdanovic 2014, 
Hayasaki et al. 2015, Shiokawa et al. 2015, Piran et al. 2015, Bon- 
nerot et al. 2015). However, this explanation begs the question: if 
circularization is so inefflcient as to result in such a massive and 
long-lived reprocessing layer, then why do observed TDE light 
curves appear to so faithfully track the predicted oc decay in 
the fail-back rate? Miller (2015) proposed that the low tempera¬ 
tures of TDE flares result from line-driven winds from an extended 
accretion disk, which suppresses the fraction of the accreting gas 
which reaches the hot inner disk. However, this requires a Keple- 
rian disk that extends to thousands of gravitational radii, which the 
stellar debris in a TDE does not have enough angular momentum 
to create. Nevertheless, the idea that the fraction of the infalling gas 
that ultimately reaches small radii to accrete is reduced may well 
be correct. 

A deflning feature of TDE fall-back is the very low angular 
momentum of the infalling gas compared to its gravitational bind¬ 
ing energy (Coughlin & Begelman 2014). Gas inflowing at super- 
Eddington rates cannot radiatively cool on the dynamical or accre¬ 
tion timescale, in which case photons are trapped and advected with 
the fluid. Low initial binding energy coupled with this globally adi¬ 
abatic evolution suggests that only a fraction of the infalling debris 
is likely to remain gravitationally bound to the BH, with the re¬ 
mainder unbound in an outflow (Strubbe & Quataert 2009; Lodato 
& Rossi 2011). 

This work develops an alternative model for thermal TDE 
emission, and the source of the reprocessing material, based on 
the assumption that the accreted fraction of gas is very low and 
that most of the disrupted star that begins weakly bound to the BH 
is promptly ejected soon after reaching pericenter. Such a massive 
outflow is necessarily slow and, in contrast to the flnding of past 
work (Strubbe & Quataert 2009), is found to remain opaque to hard 
EUV/soft X-ray radiation from the inner disk for an extended pe¬ 
riod of months or longer. The “reprocessing” of high frequency ra¬ 
diation from the inner disk due to selective absorption by this neu¬ 
tral outflow will regulate the effective temperatures of some TDE 
flares to be ~ few 10"^ - 10^ K, helping to explain the low and 
temporally constant colors of the optically-discovered events. Our 
results also provide insight into a potential unification of thermal 
TDE phenomenon, by delineating the conditions under which X- 
rays from the inner disk will promptly ionize their way out through 
the outflow ejecta, resulting instead in a luminous X-ray flare, at 
least for some range of viewing angles. 

This paper is organized as followed. In §2 we briefly review 


essential elements of TDE physics. In §3 we describe the proper¬ 
ties of the unbound outflow, and of the smaller fraction which is 
accreted by the BH. In §4 we present a toy model for the evolu¬ 
tion of the hard (ionizing) and soft (optical) radiation behind the 
ejecta. This model describes the physical process by which the ac¬ 
cretion power is converted to optical radiation and quantifies the 
radiative efficiency of TDEs. Section §5 describes the ionization 
state of the ejecta, exploring over what range of conditions hard- 
UV/X-ray photons are trapped, with implications for the diversity 
of thermal TDE flare phenomena. In §6 we describe the tempera¬ 
ture of the reprocessed emission. In §7 we discuss an implications 
of our model for the cosmological evolution of the mass and spin 
of supermassive BHs, and discuss how high measured spin values 
provide tentative evidence for a low accretion efficiency in TDEs. 
In §8 we discuss our results and summarize our conclusions. Eig- 
ure 1 provides a schematic illustration of our model, the pieces of 
which are developed as we proceed. 


2 TIDAL DISRUPTION EVENTS 

A star of mass - m^Mq and radius is tidally disrupted if the 
pericenter radius of its orbit, Rp, becomes less than the tidal radius 

R, = « 7 X 10‘^ cm 

( 2 ) 

where we have assumed a mass-radius relationship i?* ~ tnl^R^ 
appropriate to lower main sequence stars (Kippenhahn & Weigert 
1990). The orbital penetration factor is defined as yS = Rt/^p > 
with a maximum value allowed for disruption outside the horizon 
of 


ySmax ~ 12my‘^M;f (3) 

where Ribco is the radius of the innermost bound circular orbit. The 
maximum BH mass (ySmax = 1) resulting in a disruption is thus 

« 4 X lQ'M^m’l^\R;^,JAR^r^l\ (4) 


which depends on the BH spin a through Ribco(^?), where, e.g., 
?^ibco = 4Rg for a non-spinning BH. This approximate treatment 
of the highly relativistic regime near the horizon is in reasonable 
agreement with the relativistic hydrodynamics models of Ivanov & 
Cherny akova (2006). 

Disruption binds roughly half the star to the BH by a specific 
energy |Etl ~ GM^R^jR^ corresponding roughly to the work done 
by tidal forces over a distance ~ R^. The most tightly bound matter 
falls back to the BH on the characteristic fail-back timescale set by 
the period of an orbit with energy 

(5) 


which is generally, to first order, independent of the penetration 
factory6 (Stone et al. 2013; Guillochon & Ramirez-Ruiz 2013). The 
resulting rate of mass fail-back at times t » tfb is given by 

The peak Eddington ratio of the fail-back rate is given by 


4^fb(^fb) 

^Edd 
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where rj - 
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Figure 1. Schematic illustration of TDE inflow-outflow geometry and thermal emission. Stellar debris following the disruption returns to pericenter at a rate 
Mfb (dark blue). Due to the low gravitational binding energy of the debris and its inability to cool radiatively (eq. [10]), only a small fraction /in ~ 0.01 - 0.1 
joins the inner disk and accretes onto the black hole (red). The vast majority is instead unbound in a wind (blue) of velocity Vej ~ few 10^ - 10^ km s“^ 
(eqs. [13], [14]), which accumulates in an outer shell of radius /?ej and mass Mej ~ M^/2 (light blue). Extreme EUV/X-ray radiation from the inner accretion 
disk (redphotons) is initially trapped behind the expanding wind ejecta in a hot nebula of radius Rin (pink). Nebular radiation ionizes the inner wind to a depth 
A (dot-dashedyellow line), which remains less than the ejecta thickness for several fall-back times (tion > few tf^; eq. [44]) for massive BHs (M. < 10^M©; 
Fig. 5). This ionized layer is heated to a temperature of ~ few 10^ K, similar to other astrophysical environments with similar ionization parameters such 
as the ejecta in classical novae. Optical radiation from this layer more easily escapes from the ejecta due to the lower electron scattering (yellow photons), 
producing a luminous flare. Lower densities perpendicular to the plane of disruption may allow X-rays to escape from this nebula preferentially along these 
polar directions, resulting in a possible viewing angle dependence of the thermal X-ray luminosity. Broad emission lines of hydrogen or helium (dashed yellow 
photon) may originate from the photo-ionized region, or from larger radii in the inner wind or outer shell. 


3 INFLOWS AND OUTFLOWS 

Stellar debris on a nearly parabolic orbit returns to the BH with an 
angular momentum sufficient to circularize at a radius - 2R^ = 
2Rt/3~^. The poorly understood process by which matter dissipates 
its orbital energy is an area of active research (Kochanek 1994; 
Guillochon & Ramirez-Ruiz 2013; Hayasaki et al. 2015; Shiokawa 
et al. 2015). In this work we assume that circularization of at least 
a fraction of the debris is efficient on timescales ~ ^fb- We further 
assume that the narrow infalling tidal streams are sufficiently dense 
to reach the disk by penetrating the outflowing matter distributed 
over a larger solid angle (Fig. 1). 

We assume also that only a small fraction /n «: 1 of the 
gas that returns to pericenter becomes part of the Keplerian disk. 


Gas infalling at rates well above the Eddington rate (eq. [7]) 
cannot radiatively cool on the dynamical time (e.g., Strubbe & 
Quataert 2009). Quantitatively, the timescale for radiative cooling 
of gas inflowing with density p ^ Mfb/(4:7rVffRp) near pericen¬ 
ter can be estimated as the diffusion timescale t^ ~ Tes(Rp/c) ^ 
^fb^es%n/(47rcRp), whcrc Vff ~ Rp/Uyn is the free-fall velocity and 
Tes ~ pRp/^es is the electron scattering optical depth near pericenter. 
The ratio of the cooling timescale to the dynamical timescale, 


tdyn 
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thus exceeds unity until the “cooling” time 

4ooi/?fb ~ » 1, 


( 8 ) 


(9) 
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where we have used equations (2) and ( 6 ). Note that equation (9) is 
a lower limit on the cooling time, since in equation ( 8 ) is derived 
assuming the density of a spherical inflow and the inflowing tidal 
debris will be much denser. 

At times t t^ooi, gas returning to pericenter is effectively 
adiabatic on the dynamical timescale.^ On the other hand, the 
specific energy that must be removed to circularize the debris, 
\Ec\ - GM,/2/?p, is typically^ two orders of magnitude larger than 
the orbital binding energy of the debris \Et\ ~ GM.R^/R^, 

EJE, ~ R^p-^IR, ^ . (10) 

Hydrodynamical simulations by Li et al. (2013) explore such 
tenuously gravitationally bound (“zero energy”), radiatively ineffi¬ 
cient accretion flows in the context of low luminosity AGN. They 
And that only a fraction /in ~ cr of the inflowing matter becomes 
bound to the disk, where a ^ 0.01 - 0.1 is the Shakura Sunyaev 
effective viscosity parameter. This finding can be understood as a 
comparison between the accretion (or viscous) time for a thick disk, 
Lise ~ (^^Kcr)"^ and the time the tidal streams spend near pericen¬ 
ter Lrb ~ ’ where Gk is the Keplerian angular frequency. The 

remaining fraction 1 - /n ~ 1 is ejected in an outflow from near the 
circularization radius (see below). 

There are differences between the TDE problem and the set-up 
of Li et al. (2013), who evaluate a rotating flow that becomes spher¬ 
ical at large radii and who include no radiation transport. Neverthe¬ 
less, a small value of /n ~ 0.01 - 0.1 appears natural when one com¬ 
bines the requirements of ( 1 ) energy conservation; ( 2 ) globally adi¬ 
abatic evolution, due to inefficient radiative cooling on timescales 
short enough to explain observed TDE light curves; (3) matching 
the measured velocity line widths in TDE flares (see below). 

As an explicit example, consider an alternative picture to the 
‘prompt ejection’ scenario outlined above based on Li et al. (2013). 
Imagine that gas does first circularize into a circular rotating con¬ 
figuration without becoming unbound on the dynamical timescale. 
The very low binding energy of the circularized matter implies 
that it will form a geometrically thick torus with substantial radial 
pressure support (e.g., Coughlin & Begelman 2014). If a fraction 
T] of the accreted rest mass is released on the accretion timescale 
~ tdyn/o^, the total energy released will be sufficient to unbind the 
rest of the gas on a timescale Lut ~ where the specific 

binding energy Eb = niax[|EtL l^orbll is the maximum of that due 
the tidal disruption process, Et, and that of the disrupted star’s or¬ 
bit, Eorb- In this case, the condition for non-radiative evolution on 
the timescale ~ Lut is that the ratio 


A 

Lut 


T]ac^ td 
Eh ^dyn 




k 

kyn 


( 11 ) 


^ Comparing the cooling and dynamical times near pericenter is not the 
relevant criterion for energy released during the circularization process at 
larger radii, e.g. near stream apocenter. However, since the total energy dis¬ 
sipated in the circularization process is oc 1 /r and hence will be dominated 
by the energy released on the smallest radial scale ~ Rp, equation (8) even¬ 
tually becomes the relevant criterion. 

^ Speculatively, the value of Et could be dramatically increased for rare 
deeply penetrating events; at high/, the interaction between tidal compres¬ 
sion and misaligned BH spin may be capable of binding a large fraction 
of the debris much more tightly to the BH (Stone et al. 2013), but this has 
not been simulated numerically. Recent numerical results indicate such an 
enhancement may be possible for extreme / even absent BH spin (Evans 
et al. 2015), although such an enhancement can also arise spuriously in hy¬ 
drodynamical simulations of TDEs that under-resolve the orbital midplane 
(Guillochon et al. 2009). 


exceeds unity, where in the final equality we have assumed that 
Eh = Et, i.e. that |Etl ^ l^orbL as is usually well satisfied. Eor 
a ~ 0.1, the condition L > Lut is similar to the L > Lxp condi¬ 
tion given above in the prompt ejection scenario (eq. [ 8 ]) for char¬ 
acteristic parameters, e.g. ~ 1. Thus, according to 

either criterion, the debris cannot efficiently cool for several fall¬ 
back times. This renders outflows likely, except in the unlikely oc¬ 
currence that the majority of the excess thermal energy is advected 
across the BH horizon (Narayan & Yi 1995). 

Our assumption that a large fraction of the inflowing tidal de¬ 
bris is promptly unbound appears to be incompatible with the find¬ 
ings of recent hydrodynamic simulations of debris circularization 
in TDEs, which do not show such mass-loaded outflows despite 
their use of adiabatic equations of state (Hayasaki et al. 2015, Sh- 
iokawa et al. 2015, Bonnerot et al. 2015). These simulations do 
not, however, include magnetic fields, which may be responsible 
for the angular momentum transport (viscosity) that could prefer¬ 
entially bind some mass at the expense of unbinding the rest. Per¬ 
haps more importantly, for numerical reasons some of these works 
simulate the disruption of stars on unrealistically low eccentricity 
orbits, for which the initial gravitational binding energy of the star 
is much larger than the more physical case of a nearly parabolic 
(zero energy) orbit. This is true for all the simulations of Hayasaki 
et al. (2015), and some of Bonnerot et al. (2015), where eccen¬ 
tricity ^ ^ 0.8 and Eh/Ec ~ 0.1. This artificially high ratio renders 
these simulations at best marginally capable of depositing sufficient 
energy to drive an outflow from scales R ^ Rp^. Shiokawa et al. 
(2015) initializes a parabolic orbit, but they employ an artificially 
small BH to star mass ratio of 500:1, which also results in a large 
initial binding energy of Et/E^ ~ 0 . 1 , according to equation ( 10 ). 
Although enough energy is released in shock dissipation to in prin¬ 
ciple unbind some gas from the BH, this energy is quickly advected 
inward through an artificially large inner boundary (Shiokawa et al. 
(2015), their Eig. 14). 

Even if these hydrodynamic simulations were run for more 
physical mass ratios and initial orbital eccentricities, the issue of 
outflows might still not be satisfactorily resolved. The importance 
of unbound outflows in radiatively inefficient accretion flows is an 
issue of long-standing debate, even in cases when the disk begins 
tightly gravitationally bound near its outer edge (Stone et al. 1999, 
Blandford & Begelman 1999; Quataert & Gruzinov 2000; Hawley 
& Balbus 2002; see Yuan & Narayan 2014 for a recent review). 
Global MHD simulations (e.g., Hawley & Balbus 2002) And that 
a large fraction of the disk mass is lost to disk outflows within a 
decade of radius of its outer disk, with a velocity close to the escape 
speed from this location. TDE fail-back debris is even more weakly 
bound (lower Bernoulli parameter) than the initial conditions used 
in these simulations. 

The small fraction of Mfh which is added to the inner Keple¬ 
rian disk feeds the BH, producing an accretion luminosity 

Lncc = ///inMftc" ^ 2x 10^‘*ergs-‘ 

( 12 ) 

where /n-i = /m/O.Ol. The radiative efficiency of geometrically 
thin accretion varies from 77 ^ 0.04 - 0.42, depending on the spin 
of the BH and its orientation relative to the angular momentum 
of the accreting gas. Remarkably, a relatively high efficiency of 


^ Bonnerot et al. (2015) also simulate one TDE with e = 0.95 and a more 
physical Eb/Ec ~ 0.01, but run this simulation for much less than a disk 
viscous time, which (for a ~ 0.01) would be needed to launch the outflow. 
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Tj ^ 0.05 for a non-rotating BH was found by radiation (magneto- 
)hydrodynamical simulations of highly super-Eddington disk ac¬ 
cretion (e.g. Jiang et al. 2014; see also Ohsuga et al. 2005; Sad- 
owski & Narayan 2015; Jiao et al. 2015), when both radiative and 
kinetic luminosity are considered. This is contrary to the expecta¬ 
tions of previous analytic models (e.g. Shakura & Sunyaev 1973). 
In our model as well, rj is an agnostic parameter accounting for 
both radiation and the accretion power placed into the kinetic en¬ 
ergy of fast winds or jets from the inner disk. Such outflows will 
rethermalize after colliding with the inner layers of the bulk ejecta, 
in effect producing hard radiation-albeit indirectly. By assuming 
that the BH accretion rate directly tracks the fall-back rate, we are 
also neglecting viscous spreading of the accretion disk, which can 
become relevant on longer timescales (Cannizzo et al. 1990; Shen 
& Matzner 2014). 

The majority of the stellar mass is instead lost in an outflow 
from near pericenter at a rate of Mout = (1 “ ^ Mfb, re¬ 

sulting in a total ejecta mass of Mej ^ M^I2. The kinetic power of 
the unbound debris, MoutVgj/2, is that required to bind the smaller 
accreted fraction /in into the Keplerian disk, | - GM,/inMfb/2/?cL 
again under the assumption that the original energy of the debris 
is effectively zero (eq. [10]). This results in a minimum asymptotic 
ejecta velocity of 


1/2 


GM.U 

(1 ~ fn)Rc ) /in^l 


3l00kms-^ 


(13) 


For relatively high mass BHs, will approximately equal 
the mean velocity of the bulk of the ejecta, and so we take Vej = 
in most estimates to follow. However, at early times accretion 
energy deposited behind the ejecta is trapped on the expansion 
timescale (§4). Work done on the ejecta by the nebular pressure 
will in this case convert the accretion luminosity Lace (eq. [12]) into 
outflow kinetic energy, MoutVgj/2, with high efficiency. The matter 
ejected at the early times will thus be accelerated to a maximum 
velocity of 


max 


V^c « 1.3 X 10'" km s-‘ (14) 


which typically exceeds (eq. [13]) by a factor of a few. The 
observed line widths in the optical spectra of TDE flares correspond 
to expansion velocities ~ 3x10^ -10"^ km s~^ (see Table 7 of Arcavi 
et al. 2014 for a compilation), motivating values of /n ~ 0.01 - 0.1 
and 7/ ~ 0.1, consistent with the fiducial values adopted above. 

The outer radius of the ejecta grows with time as 

= v^t = 1.1 X 10*5 (15) 

where we have used equation (13). Note that exceeds by several 
orders of magnitude the outflow launching radius (eq. [2]) on the 
fail-back time. On timescales t ~ tfb over which most of the mass 
and energy are released, the outflow has plenty of time to spread 
laterally and encase the BH. 

The geometry of the wind ejecta will be complex, requiring 
more detailed numerical work to properly characterize in detail. 
However, for purposes of a simple model it is natural to decom¬ 
pose the outflow into two components: (1) a quasi-spherical “outer 
shell” comprised of most of the total mass, released mainly over 
the first few fall-back times; and (2) a quasi-steady-state “inner 
wind” comprised of matter unbound within the most recent decade 
of time since disruption. We assume that the inner wind becomes 
quasi-spherical and encases the central EUV/X-ray source by a ra¬ 
dius /?in ~ 3/?c ~ the precise value of which is also uncertain 
and will require future simulations to more accurately determine. 


A spherical geometry may not accurately characterize the disk out¬ 
flow if it preferentially emerges in one direction, such as might be 
the case if the returning debris becomes unbound during the first 
pericenter passage before “sampling” the full azimuthal extent of 
the inner disk. In §8 we address how the physical picture changes 
if the outflow is confined to a solid angle 4;7r. 

The outer shell can at times t > few tfb be approximated as a 
homologously expanding shell of mass ^ M^/2, radius Rq^, thick¬ 
ness ~ and mean density 


Psh ' 




(16) 


where we have used equation (15) for R^y The electron scattering 
optical depth is 


Tes.sh ^ Psh/^es-Rej (17) 


The shell will remain optically thick (Tes,sh >1) until a time 


tthin / ^fb 






(18) 


The density profile of the steady-state inner wind is given by 


Pw 


Mfb(l-7in) 

4;rr2vej 


3.0 X 10-“ g 

and the electron scattering optical depth is given by 



^es,w 



4^fb^es 

4;r/?inVej 




in,6 



( 20 ) 


where /?in,6 = and we have approximated the upper limit 

of integration as ~ Vejt » Rm since the value of Tes,w is not sensitive 
to this choice. 

Note that because Tes,w > ^es,sh. the outer shell will become op¬ 
tically thin to Thomson scattering earlier than the inner wind. The 
outer shell contains most of the mass, but the inner wind dominates 
the total column. This will remain true until the wind shuts off, as 
may occur once radiative cooling of gas streams inflowing to the 
disk becomes effective at t > Looi (eq. [9]) and inflowing matter 
settles into the central disk with higher efficiency. 


4 OPTICAL LIGHT CURVE MODEL 

A complete model for the coupled evolution of the nebula and the 
ejecta would require non-LTE, multi-frequency, multi-dimensional 
radiation hydrodynamics simulation, which is well beyond the 
scope of this paper. Here we instead provide a 3-zone toy model 
that captures the basic physical processes and provides a model for 
the optical light curve. 

The energy in radiation greatly dominates that of matter in 
both the ejecta and nebula. We divide this total radiation energy, 
L, into a “hard” component, Lh, and a “soft” (optical) compo¬ 
nent, Lopt. Hard radiation corresponds to EUV/X-ray photons pro¬ 
duced by the inner BH accretion disk at a characteristic temperature 
Leff ~ 10^ - 10^ K (eq. [1]), forming the nebular region around the 
BH. We assume that the soft component is produced exclusively 
by reprocessing of hard radiation after being absorbed by the ejecta 
walls. In §5 we show that at early times the ejecta is sufficiently 
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neutral to be opaque to the hard radiation, sueh that the hard ra¬ 
diation is absorbed in a thin layer separating the nebula from the 
bulk of the ejecta. This ionized layer is heated to temperatures ~ 
few 10"^ K, the radiation from which sources the optical emission 
(§ 6 ). It turns out that the inner wind is more challenging to ionize 
(requires a higher ionization parameter) than the outer shell (§5). 
Hard radiation is thus either confined to the nebula at radii > 
or it penetrates the entire ejecta and escapes from the system. 

It is assumed in this section that hard radiation is trapped be¬ 
hind the ejecta, in which case its energy evolves as 


the photon luminosity reaching the outer ejecta prior to the time 

tjtf, « 2.6^3/5^-1/2^1/5^-W^ (26) 

Combining the above, equation (23) in effect reads 


dE, 


opt 


dt 


' Aad-^ac 


-^opt 


(27) 


where the final term. 


^opt/ idi 


(28) 


^-T ^h(I-A) 

dt tic 


( 21 ) 


where Lace ^ /in^fb (eq. [ 12 ]) accounts for photon or kinetic lumi¬ 
nosity from the hot inner disk. The second term accounts for losses 
due to the absorption of hard radiation by the neutral ejecta, where 
fie = Rnic is the light crossing time of the nebula (the time between 
photon-ejecta interactions) of radius < R^^ and (1 - A) is the 
probability of absorption per interaction, where A is a frequency- 
averaged albedo. Metzger et al. (2014) show using Monte Carlo 
results (their Fig. 3) that A ~ 0.4 - 0.8 for ratios of the absorptive to 
scattering opacity, ^ ~ 0.01 - 0 . 1 , achieved in the ionized layer (see 
eq. [41] below). Hard radiation is therefore absorbed on at most a 
few light crossing times, well before adiabatic losses become im¬ 
portant. The steady state solution to equation (21), 


^H^Taccfic/(l-A), (22) 


is thus well satisfied at all times. Equation (22) will be useful later 
to estimate under what conditions the ejecta is photo-ionized by the 
nebular radiation. 

The energy in the the soft radiation evolves according to 


^^opt _ £'h( 1—A) /JLopt\ 

where the first term accounts for the absorbed hard radiation re¬ 
radiated as optical/near-UV radiation (the sink term in eq. [21]) 
and the second term accounts for adiabatic losses. Adiabatic losses 
can result from both the inner wind and the outer shell. The outer 
ejecta shell loses energy due to PdV work at a rate (JLopt/<i0ad = 
-(Lopt//?ej)Vej due to the overall expansion of the system. 

The adiabatic loss term for the inner wind is more compli¬ 
cated. Optical radiation injected at the inner edge of the wind 
(r ~ /?in) is advected outwards (“trapped”) until it reaches a crit¬ 
ical trapping radius R[^ - Mfb(l - /in)^es/( 47 rc), defined as the 
point exterior to which the photon diffusion time is shorter than 
the outflow expansion time, i.e. td/Lxp ~ Tes,wVej/c < 1 , where 
'^es.w ~ ~ Afft,/Ces/(4;rrVej) is the Thomson optical depth 

in the wind exterior to radius r (eq. [ 20 ]). 

Radiation advected from r - R^^ to Rtr has its temperature 
decreased adiabatically as L oc luminosity that 

escapes is suppressed from its injection value Lace by a factor oc 
oc oc (Strubbe & Quataert 2009). In other words, 

one can define a radiative efficiency 


r _ 

/rad = y 

t^acc 


^in 



(24) 


where L^ is the luminosity escaping from the inner wind to the 
outer shell and 

.6 ' 


\4;rcRin/ 


\ 10/9 


^fb / 


(25) 


Radiation trapping by the inner wind (/rad < 1) will thus suppress 


accounts for radiative diffusion through the outer shell and repre¬ 
sents the observed optical luminosity, where t^ ^ Tes,sh(^ej/<3) is 
the photon diffusion timescale through the outer shell. No corre¬ 
sponding diffusion term exists in equation ( 21 ) for the hard radia¬ 
tion since the opacity to hard radiation is effectively infinite until 
the ejecta becomes sufficiently photo-ionized (§5). 

Finally, the mean ejecta velocity of the outer shell evolves ac¬ 
cording to momentum conservation"^ 

^(MejVej) = (l-/i„)MftV™, (29) 

where Mej = J^'(l - /i„)Mfb(r')*' is the ejecta mass at time t. The 
source term 


V„=/radVe7 + (l-/rad)v™’^ (30) 

is the instantaneous velocity from the inner wind, which is 
smoothly interpolated between the earliest phases when the inner 
wind traps the radiation /^d 1 and Vw ~ (eq. [14]) and 

accretion power primarily goes into accelerating the outflow, and 
later phases when the inner ejecta is diffusive to radiation (/^d ~ 1 ) 
and the outflow retains its initial kinetic energy Vw ~ (eq. [13]). 
The mean asymptotic velocity of the ejecta will always lie between 
and (eq. [14]). 


4.1 Results 

Figure 2 shows an example solution of equations (27) and (29) 
for fiducial parameters: stellar mass = IM©, penetration fac¬ 
tor yS = 1, BH mass M. = lO^M©, accreted fraction /n = 0.03, 
accretion radiative efficiency 77 = 0.1, and ejecta albedo A = 0.5. 
At early times the radiated luminosity (brown line) is less than the 
accretion power of the BH (green line) as the result of adiabatic 
losses in the inner wind (/ad < 1; eq. [24]). Within a few fall¬ 
back times, however, Lrad approaches Lace to high accuracy. As will 
be discussed below, this is because once radiation escapes the inner 
wind, adiabatic losses from the outer shell are negligible. The mean 
velocity of the shell starts at Vej « > 10 "^ km s ^ because radi¬ 

ation is initially trapped, but then decreases because matter ejected 
later, after radiation freely escapes from the inner wind, has a lower 
velocity. The energy in optical radiation is larger than that in hard 
radiation by a factor ~ Tes/(1 - A) » 1. Hard radiation is absorbed 
within ~ (1 -A)“^ light crossing times, while soft radiation escapes 
on the longer diffusion time of t^ = Tesfic- 

Figure 3 shows optical light curves for the same fiducial pa¬ 
rameters as used for the solution in Fig. 2, except for a range of 
different BH masses, M, = 10^ - 10^M©. For / = 1, the optical 
flare typically lasts for several months with a peak luminosity of 


^ We neglect the BH gravitational force, which is unimportant on radial 
scales Lej » Rc of interest at times t > tfb • 
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Lpk > 10'*'* erg s"', relatively independent of M.. Figure 4 is oth- 
erwise similar, but instead shows the case of a star disrupted with 
the largest allowable penetration factor,for each BH 
mass (eq. [3]). For a deeply penetrating disruption the emission is 
seen to be strongly suppressed, despite the fact that the peak fall¬ 
back rate Mfb(?fb) oc M~^-^ is greater for small BHs, independently 
of p (eq. [ 6 ]). The light curve also peaks near the fall-back time, tfb, 
for all values of M, and p. 

The light curves in Figs. 3 and 4 can be readily understood 
from basic considerations. To aid in the following discussion. Fig¬ 
ure 5 shows key timescales as a function of M,, with other param¬ 
eters fixed at their fiducial values of 77 = 0.1 and /in = 0.03. 

We begin by addressing what causes the light curves to devi¬ 
ate from the expected oc decay at early times. First consider 
the possibility that this suppression is due to adiabatic losses from 
the outer shell. Optical radiation escapes the outer shell without 
significant losses once the radiative loss term Lj-ad in equation (27) 
exceeds the adiabatic loss term, i.e. once the ratio of the diffusion 
to the expansion timescale, 

tilt « Tes,sh(Vej/c) « Q.lmf !ta,)-\ (31) 

becomes less than unity. This condition < t is satisfied after the 
time 4 d given by 


U/tn, ^ (32) 

where we have assumed Vgj = (eq. [14]) because radiation is 
generally trapped by the inner wind on this timescale ( 4 d ^ ^tr; 
Fig. 5). The fact that 4 d < tfb for M, > lO^M© implies that adi¬ 
abatic losses from radiation passing through the outer shell does 
not suppress the peak optical luminosity for most BH masses of 
interest. 

What about adiabatic losses from the inner wind? As dis¬ 
cussed earlier, the inner wind suppresses the escape of radiation 
for a time ttr (eq. [26]). This trapping time is less than the fall-back 
time tfb (~ timescale of peak luminosity) for BH masses above a 
critical value of 


M.,„ = 6.6 X (33) 

If M. > M, tr then the peak optical luminosity is simply the 
maximum accretion power, Lpk ~ Lacc(^fb)- On the other hand, if 
M, < M, tr then ?tr > tfb and the peak luminosity will be reduced 
to Lpk = /rad^acclfpk ^ud the light curve will decay after the peak 
time as Lj-ad ^ from t = tfb until t = 4 . At times t > 4 adi¬ 
abatic losses from the inner wind become negligible, causing the 
light curve to again decay oc Lace ^ This behavior matches the 
M, < M, tr light curves shown in Figs. 3 and 4. 

Quantitatively, the peak optical luminosity is given by 

6xl0«ergs-‘yS-2/3;,_i^„ _2M^fTO0/8i?2^3^ M. < M.,„ 
~ I 2 X 10'*'* erg s"‘ ri.ifn,- 2 Mffrnf, M. > M.,,r 

where we have used equation (12) for Lace- For = 1 the peak lu¬ 
minosity is maximized for BH masses of M, = M, tr ^ 6 X 10^ Mo 
(eq. [33]; although note the sensitive dependence on the uncertain 
value of Rin) and Lpk depends weakly on BH mass for M, < M, tr 
(Fig. 3). However, for deeply penetrating TDEs with p ~ /^max 

events, the peak luminosity V for M. < M.,tr is 

more highly suppressed and reaches a maximum for higher BH 
masses (Fig. 4). 

Although the light curve peaks at t ~ tfb in all cases, the bolo- 
metric output is dominated by energy released on the maximum of 


m.= 1, /?=1, M=10®Mo, fi„=0.03, r] = 0A 



Figure 2. Model solution for TDE outflows and optical radiation as a func¬ 
tion of time since the fail-back time of the mostly tightly bound debris, cal¬ 
culated for a model with BH mass M, = lO^M©, stellar mass = IM©, 
penetration factor / = 1, BH accreted fraction /n = 0.03, ejecta albedo 
A = 0.5, inner wind radius Lin = 6 Lp, and inner disk radiative efficiency 
77 = 0.1. Quantities shown include the accretion power input from the BH, 
Lace (green; eq. [ 12 ]), the radiated optical luminosity, 4 ad (brown), mean 
velocity of outer shell, Vej (red), outer shell radius Lej (blue), energy of 
EUV/X-ray radiation in the nebula, Eu (black), and the energy of optical 
radiation in the nebula, Lopt (gold). 


tfb and 4 . The total radiated energy is thus given by 
£rad f r M. > M.,„ 

where 

^rad,max “ ~ 10 Crg 77_i/n,-2^*- (^6) 

is the maximum radiated energy. Values of Lj-ad ~ 10^^ - 10^^ ergs 
for 777 ^ ~ 1 , 77 ~ 0 . 1 , and /n ~ 0.01 - 0.1 are consistent with the 
bolometric outputs of most optical TDE flares (Gezari et al. 2008; 
Gezari et al. 2009; Gezari et al. 2012; Chornock et al. 2014; see 
Eig. 2 of Stone & Metzger 2014 for a compilation). 


5 IONIZATION STATE 


In this section we address under what conditions the ejecta will be 
opaque to hard nebular radiation above the ionization thresholds of 
hydrogen and helium. As the density of the ejecta decreases, ioniz¬ 
ing radiation from the central source will eventually penetrate this 
material and escape to infinity. After this point, hard radiation will 
be thermalized much less effectively. We consider the ionization of 
(34)both the inner wind and of the outer ejecta shell, since it is not clear 
a priori which is the most effective at trapping the hard radiation. 

Absent an external source of photo-ionization, atoms will re¬ 
combine into their ground states as the ejecta expands and adiabat- 
ically cools. Eor hydrogen-like species of charge Z, atomic mass 
A = 2Z and mass fraction A 4 , the bound-free opacity of neutral gas 
is approximately given by (Osterbrock & Eerland 2006) 


^bf 


^A^hf,vfn 

IZnip 


1.8 X 10 ^ cm^ g-‘ 



(37) 


where crbf,v is the bound-free cross section for photons of frequency 
y, Vthr ~ 13.6Z^ eV is the ionization threshold frequency and /n is 


(35) 
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Figure 3. Optical luminosity versus time since disruption, calculated for 
fiducial parameters /in = 0.03, = 1,/? = 1,77 = 0.1 and shown for 

different BH masses M./M© = 10^ {solid blue), 3 x 10^ {dotted red), 10^ 
{dashed black), 3 X 10^ {dot-dashed green), \Qp {triple-dot dashed brown). 
Note the weak dependence of the peak luminosity on BH mass for M, < 
^•,tr (eq. [33]), as given by equation (34). 



100 200 300 400 

t (days) 

Figure 4. Same as Fig. 3, but for /5 = ySmax (eq- [3]), i.e. for stars which 
are tidally disrupted on deeply plunging orbits. Note the suppression of the 
peak luminosity for low mass BHs, due to adiabatic losses from the inner 
wind. 


the neutral fraction. The value of ku above the threshold is orders of 
magnitude higher than the electron scattering opacity, to which the 
ejecta already remains opaque for months to a year (eq. [17]). Hard 
radiation therefore cannot escape to the observer unless the neutral 
fraction is reduced to a value /n 1 by nebular photo-ionization.^ 
The importance of photo-ionization is commonly quantified 


^ Bound-free opacity was neglected in previous works on TDE outflow 
(e.g., Strubbe & Quataert 2009) on the grounds that Kramers opacity is sub¬ 
dominant compared to electron scattering opacity. Kramers opacity, how¬ 
ever, applies only if the atomic level populations are in collisional ioniza¬ 
tion equilibrium and does not account for photo-ionization (see Fig. 18 of 
Ferland et al. 2013). 


Figure 5. Key timescales in days as a function of BH mass. Black solid 
line: fail-back time, ^fb(eq. [5]), over which most mass and energy is re¬ 
leased. Gold solid line: radiative timescale fcooi (eq. [9]), prior to which 
matter returning to pericenter cannot cool on the dynamical time and hence 
powerful outflows from the disk are possible. Blue solid line: timescale, 
(eq. [32]), after which radiation escapes from the outer ejecta shell without 
suffering adiabatic losses. Red lines: wind trapping time ftr (eq. [25]) before 
which optical emission is suppressed by adiabatic losses by the inner wind, 
shown separately for penetration factor = I (solid line) and [i = ySmax 
(dashed line) Green solid line: approximate time hon (eq. [44]) after which 
EUV/X-rays from the inner accretion disk ionize and escape through the 
entire ejecta for typical viewing angles. Brown lines: time rthin when the 
outer shell becomes optically thin to Thomson scattering (eq. [18]), shown 
separately for penetration factor /3 = I (solid line) and = /^max (dashed 
line) 


by the ionization parameter, which for the TDE ejecta is given by 

^ _ 4;rFH _ TtcnipE^ ^ nrUpR^L^cc 

TZej n^jVn n^jVn{l-A) 

1.5 X 10^ 

4x10^ erg cm s”' , w 

where Fh ~ Ey^clAV^ is the ionizing flux of hard radiation, = 
47tRII3 is volume of the nebula, and Fn is the nebula radius, which 
we take to equal R^j in the case of the outer shell (“sh”) (eq. [15]) 
and r > Fin for the inner wind (“w”). In the final equality we esti¬ 
mate energy of ionizing radiation using equation (22). The electron 
density ^ Pojlrrip is calculated using equation (16) in the case 
of the shell and equation (19) in the case of the wind. For typical 
parameters we see that ^ ~ 10^ erg cm s"^ but that while ^ for 
the wind is constant in time, ^ for the outer shell decreases. Also 
note that in the wind / is to first order independent of radius, im¬ 
plying that the inner edge of the wind at r ~ Fjn acts as the limiting 
step (“gate keeper”) to ionizing the entire wind. We therefore focus 
below on the conditions to ionize the wind on radial scales r ^ Fin. 

Nebular radiation will ionize the ejecta to a depth A, set by 
the location at which the effective optical depth of a photon of fre¬ 
quency y > Vthr to absorption reaches unity, i.e., 

X A 

PejKbf.v [1 + Pej'fes's] ds a Tabs(l + T^s), (39) 

where a is the depth through the ionized layer. Tabs = ^bf,yPej A is the 
optical depth through the layer to true absorption. The factor 1 -b Tes 
accounts for the additional path-length traversed by the photon due 
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to electron scattering, where Tes = Pej^esA is the electron scattering 
optical depth through the ionized layer. 

The neutral fraction /„ is determined by the balance between 
photo-ionization and recombination: 


hard radiation to escape, at a time 


tfb 


2.1X^'3Z0'5(3 - 

4.4X5'3Z0'5(3 - m0/7i; 


-3/4 

2 . 


sh 


fn - 1 + 


An r 

^rec^e J 


p^o-bf,ydy 
hr 


fn«l 


orieVn I r 


-—(Tbf,ydy 

hv 


Q^rec - y) y„/jychr(l-A)(3-r) 


^Vthr^H,Vthr^thr<^ 


^ Vthr ^thrL^QcRiitTlp 


3 X 10-976 (i--^^)(3-y) -2 


-1/3 7/30 


m. 


7x 




^ -e,4 H 7in,-2 7-l--.,6 V ?fb 

where Jy - cEyi^ylAnVy^ is the mean specific intensity of ionizing ra¬ 
diation, £^h,v ^ is the nebula radiation field per unit frequency, 
Qfj-ec = 2.6 X cm^ s~^ is the approximate case B radia¬ 

tive recombination rate, and 7^ 4 is the electron temperature in the 
ionized layer in units of 10"^ K. We define 6thr = ^H,ythr^thr/^H < 1 
as the fraction of the total ionizing radiation which resides at fre¬ 
quencies near the ionization threshold, and we have again estimated 
Eh using equation (22). We expect that 6thr ~ 1 because the energy 
of hard radiation from the accretion disk (eq. [1]) is near the He 
ionization threshold of 54 eV. 

Using equations (37) and (40), the ratio of bound-free to scat¬ 
tering opacity at the ionization threshold is given by 

^ACrbf,Vthr.^ 


^thr = — 


^bf,v 


Amp/Ces 


0.013X^Z^ (l-A)(3-y) y-0.8fi-l/2 

3 X lO-^Z^X, 


io-1/2 (—V 

%r e>4 p 4 in-2 7-1 (?fb / 

(f-A)(3-y) ^-0.8 0-1/2 I yi/-l/3 7/30 

%r p /in-2 7-6 ’ 


2/3 


Photons of frequency Vthr penetrate to a depth A which can esti¬ 
mated from equation (39) to within a factor of 2 as (Metzger et al. 
2014) 


A//?ej - (2Tes)' 


-1 ^1 _i_ /I r-1 


I+44J-I 


^thr«l 


- 1 / 2-1 


(42) 


using the value of ^thr calculated for the unattenuated nebular spec¬ 
trum (eq. [41]), where /?ej here represents either the thickness of the 
outer shell or the characteristic radial scale ~ Ri^ of the inner wind. 

For the disruption of a main sequence star, the ejecta is com¬ 
posed primarily of hydrogen (Z = 1, A = 1, ~ 0.7) and helium 

(Z = 2, A = 4, Xa ~ 0.5), in which case equation (41) shows that 
^thr < 1 for typical parameters, e.g. 6x ~ 1 and A ^ 0.5 (see Fig. .3 
of Metzger et al. 2014). Equation (42) then shows that the ejecta 
will thus remain opaque to ionizing radiation (A < /?ej) until the 
electron scattering optical depth decreases below a critical value^ 
of 


, /-1/2 

' ^thr 

19- 


p/2 


,_pl/4 f. 

p/2z3/2(3 )1/27^ Jv 

1/2 

40 .., 




ii) 


-1/3 


xy^z3/2(3-7)l/2 


1/2 

2d- 


1 K'P 


where we have taken fiducial values of A = 0.5 and = 3 x 10"^ K, 
representative of the temperature in the photo-ionized layer (§6). 
The ejecta will be completely ionized by the inner disk, allowing 


2/3 


sh 

(43) 


-1 where we have used equations (17) and (20) for Tes,sh(0 and Tes,w(0» 
respectively. 

Equation (44) shows that ionizing radiation usually takes 
longer to tunnel through the inner wind than the outer ejecta shell. 
Eor fiducial parameters of M, ^ lO^M©, - I, j3 - rj = 0.1, 
1 /in = 0.03 and a spectral slope 7^1/3 expected for the multi- 
(4b)color blackbody inner accretion disk, the ejecta should thus remain 
opaque to > 13.6 eV radiation (Z = 1; X^ = 0.7; 6thr ~ 0.1) and to 
> 54 eV radiation (Z = 2; X^ = 0.5; 6thr ~ 0.5) for several months 
(Fig. 5). 

Equation (44) shows that ionizing radiation will escape at ear¬ 
lier times as compared to a case corresponding to our fiducial pa¬ 
rameters for (1) a larger accreted fraction /n; (2) higher radiative 
efficiency of the inner disk 77; (3) less massive stars < 1. Note, 
however, that the product of /n// is constrained by the observed ra¬ 
diative efficiency of the flares according to equation (35) to near 
their fiducial values. Radiation will also escape earlier than esti¬ 
mated by equation (44) along directions through the ejecta with 
significantly lower than average density, e.g., along the axis per¬ 
pendicular to the plane of the disrupted orbit. The solid angle over 
which X-rays are free to escape will increase with time, approach¬ 
ing unity at t > tion- 

Figure 6 illustrates the parameter space of possible TDE ther- 
, mal emission in the space of BH mass and penetration factor p for 
(41jfiducial parameters of 77 = 0.1 and /n = 0.03. The upper right cor¬ 
ner of the diagram corresponds to events with p > yS^ax (eq. [3]), 
for which the star is swallowed whole by the BH before being dis¬ 
rupted. For/5 < ySmax thcrc are two regimes. For M, < lO^M© ioniz¬ 
ing radiation is trapped behind the inner wind for at least a couple 
of fall-back times (dashed blue line), reprocessing the emission to 
lower frequencies and producing an optically luminous flare. Con¬ 
tours of total radiated emission are shown as dot-dashed lines. For 
M. > Mq the ejecta becomes ionized within a couple fall-back 
times (tion < 2tfb; dashed blue line; eq. [44]), in which case much 
of the bolometric output should escape at EUV/soft X-ray wave¬ 
lengths, producing an X-ray luminous flare. 

In reality, the precise division between optically and X-ray lu¬ 
minous is theoretically uncertain given our simplistic treatment of 
ionization layer physics and the likely role played by viewing an¬ 
gle. Shown for comparison with circles and X-rays near the top of 
Fig. 6 are the inferred BH masses for the current sample of ob¬ 
served TDE flares, including some events detected as both opti¬ 
cal and X-ray sources (shown with overlapping symbols). A weak 
trend of higher BH masses associated with X-ray luminous TDEs 
is apparent, although some low mass BHs also host X-ray luminous 
sources, which in our model would need to result from a fortuitous 
viewing angle. It should also be noted that large systematic uncer¬ 
tainties exist associated with BH mass measurements, which here 
are all taken from galaxy scaling relations. 


6 TEMPERATURE OF THE OPTICAL RADIATION 


^ Equation (42) provides the penetration depth of photons at the threshold 
frequency Vthr. Photons with frequencies above the threshold penetrate to a 
distance which is at most a factor of ^ 2 - 3 times larger than those at the 
threshold frequency (Metzger et al. 2014). 


Here we address the temperature of the emerging optical/UV radi¬ 
ation at times when the EUV/X-ray photons are absorbed by the 
neutral ejecta. TDE nebulae are characterized by ionization param¬ 
eters of / ~ 10^ - 10^ (eq. [38]), which are similar to those in 
other astrophysical environments. One example is the ejecta from 
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Figure 6. Thermal emission from TDEs in the parameter space of BH mass 
M. and penetration factor /?, calculated for the disruption of a solar mass 
star = M© and fiducial parameters: /in = 0.03, rj = 0.1. For deeply 
penetrating orbits (fi > ySmax; solid black line, eq. [3]), little or no emission 
is expected because the star is swallowed whole by the BH before being 
tidally disrupted (shown separately for the cases of a non-spinning BH and 
a maximally spinning BH). For BH masses of M, < M, tr (solid red line; 
eq. [33]) thermal emission is suppressed because radiation is trapped by 
the inner wind on the fail-back time. For the highest BH masses that allow 
tidal disruption. X-rays from the inner disk ionize the ejecta on timescales 
< 2ffb (dashed blue line; eq. [44]), allowing this emission to escape (“X- 
ray Luminous”). Shown with dot-dashed red lines are also the contours of 
constant bolometric radiated optical energy (eq. [35]). In addition to the 
dependence on M. and highlighted by this diagram, some events may 
be X-ray luminous even in the “Optically Luminous” parameter space for 
viewing angles with lower than average line-of-sight density, or at late times 
» 2ffb once the ejecta has been ionized. Shown with circles and crosses are 
the estimated BH masses of those TDE flares detected in optical and/or soft 
X-ray frequencies, respectively (Gezari et al. 2006, 2008, 2009; Cenko et al. 
(2012a), Arcavi et al. 2014; Chornock et al. 2014). 


a classical nova, when it is irradiated by the supersoft X-rays from 
the heated WD surface on timescales of months to years after the 
nova eruption.^ Photo-ionization under such conditions heats the 
gas to a temperature ~ 2-5x1that depends only weakly 
on the specific radiation field due to the sensitive dependence of 
the rate of line cooling on temperature (see, e.g. Cunningham et al. 
2015; their Fig. 3). The hard radiation absorbed by the ejecta should 
thus be reprocessed into emission lines and thermal radiation with a 
characteristic energies similar to the observed temperatures of TDE 
fiares. 

If optical photons produced in the ionized layer enter equilib¬ 
rium with matter as they diffused out of the ejecta, then the temper¬ 
ature of the radiation can be estimated as 


T^TE 
^ efF 


Atto-RI, 

rJ 


1/4 


: 2.1x10^ K 


t 

^fb j 

(45) 


- 11/12 


where we have used equation (15) for and in the second line 


we have assumed t > ?tr (eq. [25]), i.e. that the radiated luminosity 
tracks the accretion power, Lj-ad = ^acc ^ 

Equation (45) predicts characteristic temperatures within a 
factor of a few of those of observed optical TDE fiares. However, 
it also predicts that the temperature should decline in time, in con¬ 
tradiction with the temporally constant temperatures of optically- 
discovered TDE fiares. TDE spectra, unlike those of most super¬ 
novae, do not evolve to become redder with time as the ejecta ex¬ 
pands and cools (e.g., Chornock et al. 2014). 

The true color temperature of the emission will, however, re¬ 
main higher than that predicted by equation (45) because as the 
ejecta spreads out and becomes optically thin, the photosphere 
will recede deeper than our naive estimate of used above. 
Note in particular that for M, > lO^M©, the outer ejecta shell 
becomes optically Thomson-thin within less than a few fail-back 
times (Eig. 5). The true photosphere radius will thus increase with 
time less rapidly than oc t (and may even decrease), resulting in 
a much slower evolution of Teff- Due to the lower scattering op¬ 
tical depth of the outer layers, the mean photon energy of ~ few 
10"^ K set by the temperature of the inner photo-ionized ejecta may 
well be preserved in the emission that escapes to the observer, even 
though the ejecta is expanding as a function of time. Euture work 
is required to more accurately estimate the effective temperature of 
the reprocessed emission. 


7 IMPLICATIONS FOR COSMIC EVOLUTION OF 
SUPERMASSIVE BLACK HOLE SPIN 

Our model has significant implications for the mass and spin evo¬ 
lution of low-mass BHs. Theoretically (observationally) calculated 
TDE rates V ~ 10""^ (10"^) gal"^ yr“^ (cf. Stone & Metzger 2014, 
van Velzen et al. 2011) imply that BHs with M, < 10^M© (10^M©) 
should see significant growth from star capture over a Hubble time, 
and may even place a lower limit on BH mass, provided the tradi¬ 
tional assumption that /n ~ 1 is true. 

In our model, /n 1 and TDEs are unable to contribute to the 
growth of BHs in typical galaxies. In the traditional picture (/n = 
1), TDEs can contribute to spin evolution even at larger (~ lO^M©) 
BH masses. The Thorne limit of a. = 0.998 (Thorne 1974) is often 
taken to represent the maximum spin of BHs in the universe^, and 
BHs that grow through the accretion of coplanar gas in thin disks 
will eventually reach it (Berti & Volonteri 2008). However, because 
the set of tidally disrupted stars has little net angular momentum, 
their net effect over time is to despin rapidly rotating BHs, as shown 
in Eigure 7. If TDEs result in the accretion of a mass M^jl of gas 
onto the BH, then even a relatively low TDE rate can impose a 
different spin limit on BHs: 

«max ^ 1 - (46) 

where we have assumed isotropic final orbits for the disrupted 
stars. A more accurate estimate of ^max can be achieved using 
Eqs. 5.267 of Merritt (2013), which describe exactly the relativis¬ 
tic, inclination-dependent change in BH spin following accretion of 
stellar matter; however, the simple Newtonian estimate of equation 


^ The velocity of nova ejecta is typically a few thousand km s“^ similar to 
that of TDE outflows. The ejected mass of ~ 10“^ - 10“^M© (and hence 
density) is typically a factor of 10^ - 10^ lower than in TDEs, but the ioniz¬ 
ing luminosity of ~ 10^^ - 10^^ erg s“^ is typically smaller 10^-10^ times 
smaller, resulting in a similar ionization parameter. 


^ However, the maximum spin that can be reached through constant pro¬ 
grade accretion depends on accretion mode, and the Thorne limit strictly 
applies only to thin disks. Most optimistically, S^dowski et al. (2011) find 
that super-Eddington and low-viscosity slim disk accretion can saturate at 
<3max = 0.9995; see references therein for summaries of other models that 
find smaller saturation spins. 
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(46) is sufficient to highlight the large role TDEs can play in the 
spin evolution of smaller supermassive BHs. 

Even low, observationally inferred, TDE rates are capable of 
spinning 10^M© (10^M©) BHs down to an ^max ~ 0.8 (0.98) in the 
traditional picture of tidal disruption. Our model, however, predicts 
that the TDE contribution to BH spin evolution is vastly weaker at 
all masses, though not totally negligible. 

Strictly speaking, our assumption of isotropic final orbits for 
tidally disrupted stars is only valid for TDEs sourced from two- 
body relaxation in a spherical, non-rotating cluster; in this limit, a 
BH that grew only through the capture of N stars would statistically 
approach (a.) ^ 1/ ^[N. If the stellar cusp possesses net rotation, 
tidally disrupted stars will also, and Eq. 46 would be modified to 

^max ~ 1 -A^(l -/pin);r/in^HM^/(2M.)-A^/pin/intHM^/(2M.), where 
X is the fraction of stars on retrograde orbits and /pin is the fraction 
of TDEs from the “pinhole” regime of disruption^; because /pin ~ 1 
when M. < 10^M© (Stone & Metzger 2014), we neglect the effects 
of stellar cusp rotation^^. Our isotropic approximation is further 
justified by the fact that two-body relaxation may not be the dom¬ 
inant dynamical process for creating TDEs. TDE rates can be en¬ 
hanced above the level set by two-body relaxation if axisymmetry 
(Magorrian & Tremaine 1999; Vasiliev & Merritt 2013) or triaxial- 
ity (Merritt & Boon 2004) exist in the stellar potential; in this case, 
orbits will not conserve all components of angular momentum, and 
some will collisionlessly wander to very small pericenter. These 
aspherical channels for producing TDEs primarily funnel stars to 
the BH in the “pinhole” regime of disruption (Vasiliev & Merritt 
2013), which severely reduces the net spin they can impart to the 
BH. 

In practice, TDEs have the greatest impact on ^max if BHs 
grow through “coherent” gas accretion, i.e. accretion of gas from 
a fixed orbital plane. Whether or not this occurs depends on the 
larger-scale physics of the host galaxy. State of the art cosmological 
simulations (with subgrid models for BH spin evolution) find that 
the small BHs which likely dominate TDE rates (10^ < M,/M© < 
10^) spin up efficiently through gas accretion, generally coming 
close to a, = 0.998 (Dubois et al. 2014, Eigure 7). As shown in our 
Figure 7, this mass range is one in which TDEs render the Thome 
limit unattainable if /n ~ 1, and the true ^max at lO^M© is between 
0.6 and 0.98, depending on the TDE rate. Roughly 20 supermas¬ 
sive BH spins have been measured using the Fe Kcr technique, and 
tantalizingly, several fall into the region of parameter space that 
can constrain /n. MCG-6-30-15 has M. = 2.9^} ^ x 10^M© and 
a. > 0.98, Ark564 has M. ~ 1.1 x lO^M© and a. = 0.96+goj, 
IRAS 13224-3809 has M. ~ 6.3 x 10^ and a. > 0.987, NGC 4051 
has M. = 1.91 + 0.78 x lO^M© and a. > 0.99, and 1H0707-495 has 
M, ~ 2.3 X 10^M© and a. > 0.97 (Reynolds 2013, Table 1). 

All of these mass-spin combinations are impossible if /n = 1 
and their host galaxies have TDE rates in line with the predictions 
of Stone & Metzger (2014). Several are even more constraining; for 

^ In the pinhole regime of tidal disruption, the per-orbit change in specific 
angular momentum |A/| » |/lcI = ^|2GM^Rl, the specific angular mo¬ 
mentum of the loss cone. In this regime the star wanders in and out of the 
loss cone many times per orbit, and its final, disruptive passage through the 
loss cone is just as likely to be prograde as it is to be retrograde with respect 
to any reference axis. In the opposite, “diffusive,” regime, |A/| |/lcI and 
disrupted stars retain a memory of the cluster rotation axis. 

Beloborodov et al. (1992) examined the equilibrium spins of supermas¬ 
sive BHs that grew through tidal disruption, but assuming that the tidally 
liberated gas only accreted after settling into the cluster equatorial plane. In 
our terminology, this is equivalent to the extreme case /pin = 0 and;^ = 0. 



M. (M©) 

Figure 7. The maximum attainable supermassive BH spin, flmax, for BHs 
of a given mass M. which have grown partially through stellar tidal dis¬ 
ruption. We plot 1 - flmax to show the variation in very high values of 
BH spin. Red, blue, and black lines assume /n = 1, 0.1, and 0.01, re¬ 
spectively. The solid lines show an observationally estimated TDE rate 
of N = 10“^ gar^ yr“^ (van Velzen & Farrar 2014), while the dashed 
lines show the (higher) theoretically predicted TDE rate of A = 2.9 x 
10-5 gaj-1 yj-i from Stone & Metzger (2014). All 

curves asymptote to the Thorne limit (^max = 0.998) once the BH mass be¬ 
comes large enough that star capture is a negligible contribution to the BH’s 
spin evolution. In the traditional picture, where /m = 1, TDEs place a more 
restrictive cap on BH spin than does accretion physics for M. < lO^M©, 
but this is not the case in our revised TDE model with /m ^ 1. Shown with 
circles and arrows are published measurements and upper limits, respec¬ 
tively, on BH spin, compiled by Reynolds (2013, Table 1). We assume the 
accretion of stars with an average mass = 0.3M©. 

example, the measured mass-spin pair in NGC 4051 cannot be pro¬ 
duced by a galaxy following our theoretically predicted TDE rates 
unless /in <0.1. Although this comparison with data is suggestive, 
we caution that /n can only be rigorously constrained through sta¬ 
tistical analysis of a large sample of measured supermassive BH 
spins, a task beyond the immediate scope of this paper. 


8 DISCUSSION AND CONCLUSIONS 

We have presented a model for the thermal optical emission from 
TDEs, with implications for a possible unification of thermal TDE 
phenomenon. Our main assumption is that only a small fraction 
/in ~ 0.01 - 0.1 of the initially bound stellar debris ends up reach¬ 
ing the inner circular disk and accreting onto the BH, with the re¬ 
mainder (comprising roughly half of the original star) blown out 
in a wind or unbound stream (Strubbe & Quataert 2009; Lodato & 
Rossi 2011) that rapidly comes to encase the BH. 

Such a low accretion efficiency is suggested by numerical sim¬ 
ulations of adiabatic accretion with low angular momentum and 
negligible initial binding energy (Li et al. 2013), which find that 
only a fraction ~ cr of the infiowing gas is accreted by the central 
BH. The precise value of the effective viscosity a depends on the 
saturation level of MHD turbulence, which could well be different 
for shear in the eccentric gaseous stream of a TDE than in a circular 
accretion disk (Guillochon & Ramirez-Ruiz 2015). The efficacy of 
MHD circularization has been discounted as the primary method 
of disk circularization on the grounds that ~ l/cr » 1 orbital pe¬ 
riods are required to fully circularize the material (Guillochon & 
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Ramirez-Ruiz 2015). However, here we have emphasized that the 
energy which must be released to circularize even a small fraction 
a > 0.01 of the inflowing gas is sufficient to unbind the remaining 
fraction 1 -cr, given the weak binding energy of the debris (eq. [10]) 
and its effectively adiabatic evolution due to radiation trapping. Nu¬ 
merical simulations of super-Eddington accretion (e.g., Jiang et al. 
2014) are not directly applicable to the problem of TDE fall-back 
if they initialize a disk which is already tightly bound to the BH. 

The importance of unbound outflows in radiatively inefficient 
accretion flows remains an issue of debate (Yuan & Narayan 2014), 
so the assumption of low fin must be taken as an assumption of 
our model rather than a proven fact. However, even if only a few 
percent of the star accretes, we find that the radiative efficiency of 
reprocessing by the outflowing debris is high enough to explain the 
bolometric output of most TDE flares (eq. [35]). This challenges 
the implicit assumption of recent work that the majority of the ini¬ 
tially bound debris must indeed circularize to explain the observed 
thermal flares. Erom this viewpoint, a larger fraction of the star 
would accrete if the stellar debris were more tightly bound. This 
could occur if a sub-population of TDEs arise from very tightly 
bound stellar orbits (Hayasaki et al. 2013), or if the interaction be¬ 
tween misaligned BH spin and tidal compression in a high p event 
could dramatically increase E[ (Stone et al. 2013). However, both 
of these possibilities are speculative and unlikely to occur except in 
a small minority of TDEs. 

The small amount of energy released in binding the fraction 
fin to the disk implies relatively low velocities of a few 10^ km s“^ 
for the bulk of the material, although at early times matter can be 
accelerated by the pressure of trapped radiation to a higher veloc¬ 
ity of > 10"^ km s"^ (eq. [14]). These velocities are broadly con¬ 
sistent with the widths of the emission lines observed in the TDE 
optical spectra of- 3-10x10"^ km s"^ (Gezari et al. 2012; Arcavi 
et al. 2014). These lines may originate directly from the inner X-ray 
ionized layer of the ejecta, or from the faster outflow expected on 
larger scales. More detailed photo-ionization and radiative transfer 
calculations will be required to create more concrete predictions for 
these emission line diagnostics, and how they vary with the other 
properties of the system using considerations similar to those in 
Guillochon et al. (2014). 

One of the biggest mysteries associated with the current ob¬ 
servational sample of thermal TDE flares is what distinguishes the 
optically and X-ray luminous sources. Eor BHs of mass M, above 
a critical threshold lO^M© (which depends on fn, tj, Rin, and m*), 
we find that ionizing radiation from the inner disk is trapped behind 
the wind ejecta for timescales of several months or longer, similar 
to the characteristic duration of TDE light curves. The irradiated 
inner wind on a radial scale of > lORp can thus provide the “repro¬ 
cessing” layer needed to produce bright optical emission, without 
requiring a radially-extended distribution of mass due to inefficient 
circularization. 

Accretion power is converted into reprocessed radiation with 
high efficiency because ionizing radiation is absorbed by the neu¬ 
tral ejecta on at most a few light crossing times across the nebula. 
The diffusion timescale of the reprocessed optical radiation by the 
outer ejecta is generally short compared to the expansion timescale 
for BH masses of M« > 10^M© (eq. [32]; Eig. 5). However, for 
M. < M. tr ~ few IO^Mq (eq. [33]), trapping by the inner wind 
reduces the radiative efficiency at early times. In this case, there 
is an early phase where the radiated luminosity Lj-ad oc t be¬ 
fore the light curve approaches the the well-known oc decay 
(Eig. 3,4). The bolometric output of the TDE flare can therefore be 
considerably suppressed for low-M, and/or high p events. The tem¬ 


perature of the emission may be regulated to a relatively constant 
value of several 10"^ K by the balance between photo-ionization and 
line cooling, similar to X-ray irradiation of nova ejecta of similar 
ionization parameter (eq. [38]). 

The model developed here resembles the disk wind models of 
Strubbe & Quataert (2009) and Lodato & Rossi (2011), which have 
been discounted as the source of observed optical emission based 
on both spectral line widths and a predicted fast light curve decay 
(oc in disagreement with observations (e.g. Gezari et al. 

2012). We find that these deficiencies are remedied by consider¬ 
ing higher ejecta masses and slower outflow velocities than in past 
work (which are also more consistent with the findings of numer¬ 
ical simulations and basic energetic considerations), and by taking 
into account the role of bound-free opacity on the escape of ion¬ 
izing radiation. Our model is similar to that which has been devel¬ 
oped for optical emission from stellar explosions powered by a cen¬ 
tral engine, which is one of the most promising models for super- 
luminous supernovae (SLSNe; Kasen & Bildsten 2010; Woosley 
2010; Metzger et al. 2014). It is interesting to note that the blue 
colors and spectra of some TDE flares qualitatively resemble those 
of SLSNe (Chornock et al. 2014). 

Eor more massive BHs with M, > lO^M© (the precise value 
depending on fn and mf). X-rays from the inner disk may com¬ 
pletely ionize and escape the ejecta on a timescale tion (eq. [44]) 
exceeding a few fail-back times of the most bound debris. If and 
once this occurs, reprocessing to optical frequencies will be less ef¬ 
fective and the event will be more X-ray luminous. Such cases may 
explain the thermal X-ray detected flares (e.g.. Bade et al. 1996, 
Grupe et al. 1999, Komossa & Bade 1999, Komossa & Greiner 
1999, Donley et al. 2002, Esquej et al. 2008, Maksym et al. 2010, 
Saxton et al. 2012). 

Even in optically luminous flares where no X-ray emission is 
initially detected. X-rays may escape at late times as the density 
of the expanding shell decreases and the ionization front reaches 
the surface. We strongly encourage late-time X-ray monitoring of 
optically-discovered TDEs on timescales of months or years after 
the disruption, to search for this “ionization break-out” signal. The 
prediction that X-ray luminous TDEs may preferentially occur in 
higher mass galaxies appears to be in tension with observations of 
X-ray flares in some low mass galaxies (Eig. 6), but it is hard to 
draw firm conclusions due to both small-number statistics and the 
unreliable BH mass estimates in current TDE candidates. 

Even if X-rays are trapped by the bulk of the ejecta, es¬ 
cape is possible along directions with lower than average density 
{tion ^ eq. [44]). These viewing angles may include directions 
perpendicular to the orbital plane of the disrupted star, or oppo¬ 
site the dominant direction of the outflowing gas if the circulariz¬ 
ing debris cannot phase mix before most of the mass is unbound. 
TDEs viewed along such low-density ‘holes’ could be both X-ray 
and optically-luminous. Even if our assumption of quasi-spherical 
ejecta is not justified and ionizing radiation is not trapped behind 
the neutral outflow for more than a light crossing time, those X- 
rays that do intercept dense streams of ejecta should thermalize 
with high efficiency, given the low albedo A ~ 0.4 - 0.7 which 
is expected to characterize the exposed inner edge of the ejecta for 
the ratio of absorption to scattering opacities of relevance (compare 
eq. [41] with Eig. 3 of Metzger et al. 2014). 

Stone & Metzger (2014) highlight the tension, of at least one 
order of magnitude, between the rates of TDEs predicted by mod¬ 
eling two-body scattering of stars in galactic nuclei of ~ 10""^ gaL^ 
yr"^ (see also Wang & Merritt 2004, Merritt 2015), and the detected 
rates of optically-luminous TDEs of ~ 10"^ gaL^ yr"^ (e.g. van 
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Velzen et al. 2011). Part of this discrepancy could be the result of 
dust obseuration in galactic nuclei or observational selection ef¬ 
fects; however, part of it could be attributed to only a small fraction 
of all TDEs being optically luminous. Our finding that the peak lu¬ 
minosities and optical fiuences of TDE fiares are smaller for low 
mass BHs than would be naively anticipated from their peak accre¬ 
tion rate (eq. [34]) would aet to suppress the rate of TDE detection 
in low mass galaxies. Indeed, the inferred BH masses of optically- 
discovered TDE fiares seem to cluster around M, ~ lO^M©, despite 
the theoretical expectation that the detected TDE rate by a fiux- 
limited survey should be dominated by lower-mass BHs, assuming 
peak optical luminosity tracks the mass fail-back rate alone (Stone 
& Metzger 2014). 

Our finding that most TDEs may accrete only a small frae- 
tion of their debris makes the jetted class of TDEs all the more 
mysterious. The prototypical jetted TDE, Swift }\6AA+51, radiated 
~ 10^"^ ergs of isotropic X-rays (Bloom et al. 2011; Burrows et al. 
2011), while the kinetic energy of the jet responsible for the radio 
afterglow likewise approached a significant fraction of a solar rest 
mass ~ 10^^ - 10^^ ergs (e.g., Berger et al. 2012; Bamiol Duran 
& Piran 2013; Tchekhovskoy et al. 2014; Mimica et al. 2015). Ex¬ 
plaining this required accretion power might require a massive star 
in the eontext of low aceretion efficiency /in ~ 0.01 - 0.1, or some 
unusual circumstances, e.g. more tightly bound stellar debris due 
perhaps to an eccentric stellar trajeetory or effects of BH spin. In¬ 
terestingly, one of the consequences of our assumption of a low ae- 
creted fraction is that the true Eddington ratio of the aecretion rate 
reaching the BH is a factor of ~ 10-100 times lower than that of the 
fall-back (eq. [7]). In this case, rather low BH masses of lO^M© 
may be required to produce highly super-Eddington accretion on 
smaller scales. We speculate that if highly super-Eddington aecre¬ 
tion is the key requirement to producing a powerful jet, that this 
condition might only be satisfied by potentially rare low mass BHs. 

The potentially out-sized role of TDEs in supermassive BH 
spin evolution, which to our knowledge has not been examined be¬ 
fore, turns out to depend critically on the value of f^. For traditional 
values of /n ~ 1, it is not possible for small supermassive BHs to 
reach the theoretical upper limits on a. predicted by Thorne (1974) 
and successor works. Our prediction, that /n 1, gives TDEs 
much less leverage to cap supermassive BH spin. The Fe Kcr mea¬ 
surement of several rapidly spinning BHs with M, ~ lO^M© is hard 
to reeoncile with theoretically predicted TDE rates (Stone & Met¬ 
zger 2014) and the assumption that /n = 1. Future analyses of Fe 
Kof spin samples that account for systematic errors and small num¬ 
ber statistics should be able to rigorously constrain a combination 
of the true TDE rate and /n. 

Our model will be improved in future work in several direc¬ 
tions. First, we will employ a one or multi-dimensional outflow 
model to replace our simplified one-zone approach. Also needed 
is a self-consistent photo-ionization caleulation to more aecurately 
determine the conditions for ionizing radiation to escape, and to 
explore emission line diagnostics. 
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